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Abstract 

We build on the formulation developed in Sridhar & Singh {JFM, 664, 265, 2010), and present 
a theory of the shear dynamo problem for small magnetic and fluid Reynolds numbers, but for 
arbitrary values of the shear parameter. Specializing to the case of a mean magnetic field that is 
slowly varying in time, explicit expressions for the transport coefficients, an and ijimi, are derived. 
We prove that, when the velocity field is non helical, the transport coefficient an vanishes. We then 
consider forced, stochastic dynamics for the incompressible velocity field at low Reynolds number. 
An exact, explicit solution for the velocity field is derived, and the velocity spectrum tensor is 
calculated in terms of the Galilean-invariant forcing statistics. We consider forcing statistics that 
is non helical, isotropic and delta-correlated-in-time, and specialize to the case when the mean- 
field is a function only of the spatial coordinate X3 and time r ; this reduction is necessary for 
comparison with the numerical experiments of Brandenburg, Radler, Rheinhardt & Kapyla {ApJ, 
676, 740, 2008). Explicit expressions are derived for all four components of the magnetic diffusivity 
tensor, rjij (r) . These are used to prove that the shear-current effect cannot be responsible for 
dynamo action at small Re and Rm, but for all values of the shear parameter. 
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I. INTRODUCTION 



Astrophysical systems like planets, galaxies and clusters of galaxies possess magnetic fields 
which exhibit definite spatial ordering, in addition to a random component. The ordered 
(or "large-scale") components are thought to originate from turbulent dynamo action in 
the electrically conducting fluids in these objects. The standard model of such a process 
involves amplification of seed magnetic fields due to turbulent flows which lack mirror- 
symmetry (equivalently, which possess helicity) [l-3|. The turbulent flows generally possess 
large-scale shear, which is expected to have significant effects on transport properties jj]; 
however, it is not clear whether the turbulent flows are always helical. Recent work has 
explored the possiblity that non-helical turbulence in conjunction with background shear 
may give rise to large-scale dynamo action js-lOj. The evidence for this comes mainly from 



direct numerical simulations 



7|, but it by no means clear what physics drives such a 



dynamo. One possibility that has received some attention is the shear-current effect [10|, 
where an extra component of the mean electromotive force (EMF) is thought to result in 
the generation of the cross-shear component of the mean magnetic field from the component 



parallel to the shear flow. However, there is no agreement yet whether the sign o 
coupling is favorable to the operation of a dynamo; some analytic calculations 11 



such a 



12| and 



numerical experiments |5i] find that the sign of the shear-current term is unfavorable for 
dynamo action. 



:iow of an incom- 
who used the 



A quasilinear kinematic theory of dynamo action in a linear shear 
pressible fluid which has random velocity fluctuations was presented in ! 
"second order correlation approximation" (SOCA) in the limit of zero resistivity. Unlike 
earlier analytic work which treated shear as a small perturbation, this theory did not place 
any restriction on the strength of the shear. They arrived at an integro-differential equation 
for the evolution of the mean magnetic field and argued that the shear-current assisted dy- 
namo is essentially absent. The theory was extended to take account of non zero resistivity 
in [l^; this is again nonperturbative in the shear strength, uses SOCA, and is rigorously 
valid in the limit of small magnetic Reynolds number (Rm) but with no restriction on the 
luid Reynolds number (Re). The kinematic approach to the shear dynamo problem taken in 
uses in an essential manner the shearing coordinate transformation and the Galilean 



invariance (which is a fundamental symmetry of the problem) of the velocity fluctuations. 
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The present work extends [ij] by giving definite form to tlie statistics of the velocity field; 
specifically, the velocity field is assumed to obey the forced Navier-Stokes equation, in the 
absence of Lorentz forces. 

In section II we begin with a brief review of the salient results of [l^. The expression 
for the Galilean-invariant mean EMF is then worked out for the case of a mean magnetic 
field that is slowly varying in time. Thus the mean-field induction equation, which is an 
integro-differential equaton in the formulation of [l^ now simplifies to a partial differential 
equation. This reduction is an essential first step to the later comparision with the numerical 
experiments of jsl. Explicit expressions for the transport coefficients, an and rjimi, are derived 
in terms of the two-point velocity correlators. We then recall some results from |14|], which 
express the velocity correlators in terms of the velocity spectrum tensor. This tensorial 
quantity is real when the velocity field is non helical; we are able to prove that, in this case, 
the transport coefficient an vanishes. Section III develops the dynamics of the velocity field 
at low Re, using the Navier-Stokes equation with stochastic external forcing. An explicit 
solution for the velocity field is presented and the velocity spectrum tensor is calculated in 
terms of the Galilean-invariant forcing statistics. For non helical forcing, the velocity field is 
also non helical and the transport coefficient an vanishes, as noted above. We then specialize 
to the case when the forcing is not only non helical, but isotropic and delta-correlated-in- 
time as well. In section IV we specialize to the case when the mean-field is a function only 
of the spatial coordinate X3 and time r ; this reduction is necessary for comparision with 
the numerical experiments of jsj. Explicit expressions are derived for all four components 
of the magnetic diffusivity tensor, i]ij{T) , in terms of the velocity power spectrum; the 
late-time saturation values, rjfj , have direct bearing on the growth (or otherwise) of the 
mean magnetic field. Comparisons with earlier work — in particular jsl — are made, and the 
implications for the shear-current effect are discussed. We then conclude in section V. 

II. MEAN FIELD ELECTRODYNAMICS IN A LINEAR SHEAR FLOW 
A. Mean field induction equation for small Rm 



We begin with a brief review of the main results of |l4!]. Let (61,62,63) be the unit 
basis vectors of a Cartesian coordinate system in the laboratory frame. Using notation. 
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X = {Xi, X2, X^) for the position vector and r for time, we write the fluid velocity as 
{SXie2 + v), where 5* is the rate of shear parameter and v{X, r) is an incompressible and 
randomly fluctuating velocity field with zero mean. The mean magnetic field, B{X,t), 
obeys the following (mean-field induction) equation: 

— + SX^-—]B - SB.e^ = Vx£ + r^V^B (1) 



dr 8X2/ 

where 77 is the microscopic resistivity, and £ is the mean electromotive force (EMF), £ = 
(vxb), where v and b are the fluctuations in the velocity and magnetic fields. 

To lowest order in Rm, the evolution of the magnetic field fluctuations, now denoted by 
b^^\ is governed by, 

(|: + ^^1^) b^'^ - Sbf^e2 = iB.V)v - iv.V)B + r^V'b^'^ (2) 

This equation was solved by making a shearing-coordinate transformation to new spacetime 
coordinates and new field variables. The new spacetime variables, {x,t), are given by, 

xi = Xi , X2 = X2 - StXi , X3 = X3 , t = r , (3) 

where x may be thought of as the Lagrangian coordinates of a fluid element in the back- 
ground shear flow. The new field variables are component-wise equal to the old variables: 

H{x,t) = B{X,t), h{x,t) = 6(°)(X,r), u{x,t) = v{X,t) (4) 

In the new variables, equation ([2]) becomes. 



_ _ Sft.e, -[h-^- StH,- )n-[u.-- Stu,- j H + ,V'h (5) 

We need the particular solution (i.e. the forced solution) which vanishes at t = 0. This is 
given in component form as. 
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^ dt' J d^x' G^{X - X', t, t') [u'^i + S{t - t')5rn2u\i] X 

X [H[ - St'6i2H[] 
J* dt' j d^x' G^{x - x', t, t') [H'^i + Sit - t')5^2H[i] X 

X [u'l - St'5l2U]\ 



(6) 



The primes in H[ and mean that these functions are evaluated at {x',t'). The quantities 
H'^i and u'^i are shorthand for {dH'^/dx[) and jdu',^ I d x\) ^ respectively. Here Gn{r,t,t') is 
the resistive Green's function for the linear shear flow IJ, ll5|, which takes the form of a 
sheared heat kernel. The one property we will use is that Gn{r,t,t') is an even function of 
r. Otherwise, its spatial Fourier transform, defined by 



Gnik, t, t') = J dV Gr,{r, t, t') exp [-i k- r] 



exp 



-7] ( P(t - t') -Sh k2{t' - + Y - f'] 



(7) 



is more useful for our purposes. 

The mean EMF is given by £ = (^vxb^'^^^ = (uxh), where equation (jH]) for h should be 
substituted. The averaging, ( ), acts only on the velocity variables but not the mean field; 
i.e. (uuH) = (uu) H etc. The uu velocity correlators can be rewritten in terms of the 
vv velocity correlators; this is a useful step because the latter are referred to the laboratory 
frame. The velocity correlators have a very important property called Galilean invariance, 
which is shared by comoving observers, who translate uniformly with the background shear 
flow. If a comoving observer is at position ^ = {^1,^2,^3) at the initial time, then at a later 
time t, her location is given by. 



(8) 



Velocity fluctuations are defined to be Galilean-invariant, if and only if the statistical prop- 
erties of the fluctuations, as seen by any comoving observer, are identical to the statistical 
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properties seen in the laboratory frame; it follows that all comoving observers see the same 
statistics. There are two basic Galilean-invariant two-point velocity correlation functions, 
Qjmi and Rjm, which are defined as: 



RjUr,t,t') = ^t;,(Xe(^,t),t)t;™(x, (9) 

Then the mean EMF is a functional of the mean magnetic field, Hi, and its first spatial 
derivative. Him = {dHi/dxm)'- 



Siix, t) = eijm dt' j dV G^(r, t, t') Cj^iir, t, t')Hi{x - r, t') - 

- j dt' J dV Gr,{r, t, t') [tiji + S{t - t')5iitij2] D,m{r, t, t')Him{x - r, t') . 

(10) 

where Cjmi and Dj^n are two-point velocity correlators, which are derived from the more 
basic two-point velocity correlators, Qjmi and Rjmi of equations 



Cjmi{r,t,t') = Q,mi{r,t,t') + S{t-t')6m2Qjii{r,t,t'), 
Djm{r,t,t') = Rjm{r,t,t') - St'6m2 R,i{r,t,t') . (11) 
Then the time evolution of the mean magnetic field is given in the new variables by. 



—- - SHie2 = Vx£ + r]V^H, 
at 

Equations ( IT^ and (ITU]) form a closed system of integro-differential equations, determining 
the time evolution of the mean magnetic field, H{x,t). 
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B. The mean EMF for a slowly varying magnetic field 

The mean EMF given in equation fllUI) is a functional of Hi and Him- When the mean-field 
is slowly varying compared to velocity correlation times, we expect to be able to approximate 
£ as a function of Hi and Him- In this case, the mean-field induction equation would reduce 
to a set of coupled partial differential equations, instead of the more formidable set of coupled 
integro-differential equations given by f[T2]) and flTU]) . Sheared coordinates are essential for 
the calculations, but physical interpretation is simplest in the laboratory frame; hence we 
derive an expression for the mean EMF in terms of B{X, r). 

The first step involves a Taylor expansion of the quantitites. Hi and Him, occuring in 
equation (|T0|) for the mean EMF. Neglecting spacetime derivatives higher than the first 
order ones, we have 



Hi{x-r,t') = Hi{x,t) - r,Hi,{x,t) - (t - ^^^^ + 



Himix-r,t') = Himix,t) - (t-O ^^'^^'^^ + ... (13) 

We now use the mean-field induction equation f|T2l) . to express {dH /dt) in terms of spatial 
derivatives. Let L be the spatial scale over which the mean-field varies. When the mean- 
field varies slowly, L is large and the contributions from both the resistive term and the 
mean EMF are small, as is shown below. Let ^ and v^^^s be the spatial scale and root-mean- 
squared amplitude of the velocity fiuctuations. The resistive term makes a contribution of 
order (£/L)^Rm~^, which we now assume is much less than unity. Using equation ffTOj) . we 
can verify that Vx5 contributes terms of five different orders; {i/L)^ (£/L)(S'£/frms), {(^l 
{£/L)\Si/v,^,) and {e/L)\Se/v,^,f- These are all small when and {i/L){Se/v,^) 

are both much smaller than unity. That we must have (i/L) ^ 1 is natural from the 
familiar case of zero shear. The presence of shear introduces an additional requirement that 
{i/L){Si/vrms) ^ 1- We now define the small parameter, yU ^ 1, to be equal to the largest 
of the three small quantities, {£/LfRm'^ < 1, (i/L) < 1 and {£/ L) {Si/ Vrms) < 1- Then, 



dt 

and equations (IT^ give. 



' - S6i,H, + O(^) (14) 
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Hi{x-r,t') = Hi{x,t) - rpHip{x,t) - S{t-t')5i2Hi + 0(/x) 

HiUx-r,t') = Him{x,t) - S{t-t')5i2H^^ + 0(/^) (15) 
We substitute equation f[T^ in flTUl) to get, 

Si{x,t) = ei,mHi{x,t) dt' j d\Gr,{r,t,t') [Cjmi{r,t,t') - S{t - t') 6nCjm2{r,t,t')] 

- eijm Hip{x, t) dt' j d?r G^[r, t, t') Cjmi{r, t, t') 

- eiji Him{x, t) ^ dt' j dV G.lr, t, t') Djm{r, t, t') + ©(/i^) (16) 

The final step is to rewrite the above expression in terms of the original magnetic field 
variable, using, 

Hi{x,t) = Bi{X,t) 
-r-r , X dHi(xA) f d ^ ^ (9 \ , 

Therefore, for a slowly varying magnetic field, the mean EMF is given by, 

8, = au{T)Bi{X,T) - Vrmiir) ^^'}^'^^ (18) 
where the transport coefficients are given by, 

(yuir) = ^ijm dr' j dVG^(r,r,r') [Cjmi{r,T,T') - S{t - r') 6iiCjm2{r,T,T')] 

dr' j dV [r^ + Sr^^sri] G^(r, r, r') Cjpi{r, r, r) + 



Vimi[r) = e 



+ e,ji ^ dr' J dV G^{r, r, r') [Dj^{r, r, r') + 5r5^2/^,i(r, r, r')] (19) 

Then the mean-field induction equation ([1]), together with equations f|T8|) and f lT9|) . is a 
closed partial differential equation (which is first order in temporal and second order in 
spatial derivatives). 
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C. Velocity correlators expressed in terms of the velocity spectrum tensor 



The Galilean invariance of the two-point velocity correlators can be stated most com- 
pactly in Fourier-space. Let v{K,t) be the spatial Fourier transform of v{X,t), defined 
by 

v{K,t) = j d^Xv{X,T) exp[-iK- X] ; [K-v{K,t)]=0 (20) 

New Fourier variables are defined by the Fourier-space shearing transformation, 

ki = Ki + StK2, k2 = K2, = K3, t = T (21) 

It is proved in that a Galilean-invariant Fourier-space two-point velocity correlator 
must be of the form 

{vjiK,T)d*^{K',T')) = {27rfS{k-k')U,Uk,t,t') (22) 
where Hjm is the velocity spectrum tensor, which must possess the following properties: 

Uij{k,t,t') = U*ji-k,t,t') = U,i{-k,t',t) 

KiUij{k,t,t') = {h- St6ak2)Uij{k,t,t') = 

K!jllij{k,t,t') = {kj - St' Sjik2)Ilij{k,t,t') = (23) 
Now, the various two-point velocity correlators can be written as: 

Rjm{r,t,t') = J d^kUjmik,t,t') exp[ik-r] 

Qjmi{r,t,t') = -i y d'A; [k - St'6nk2] UjUk,t,t') exp[ik- r] 

Djmir, t, t') = j d?k [n,„(fc, t, t') - St'5^2 nji(fc, t, t')] exp [i k- r] 

Cjrrd{r, t, t') = -i y [h - St'Siik2] [Ujmik, t, t') + S{t - t')5„a^ji{k, t, t')] exp [i k- r] 

(24) 



Using the above expressions for Dj^ and Cjmi in equations (fT^ . the transport coefficients 
aui^r) and rjimiir) can also be written in terms of the velocity spectrum tensor. 
The correlation helicity may be defined as, 



Hcor{t,t') = eji^{vj{0,t)v^i{0,t')) = i / d^k[ki-St'Snk2]eij^Uj^{k,t,t') (25) 



From the first of equations fl23|) . it is clear that the real part of Iljm{k, t, t') is an even function 
of fc, whereas the imaginary part is an odd function of k. Hence only the imaginary part 
of Tljra{k,t,t') contributes to the correlation helicity. We shall see that the forced velocity 
fields we deal with later in this article possess a real velocity spectrum, and their correlation 
helicity vanishes. In this case. 



Cj^i{r, t, t') = / d=^A; [ki - St'Snk2] [n,„,(fc, t, t') + S(t - t')<5„2n,i(fc, t, t')] sin [k- r] 



are both odd functions of r. Since the resistive Green's function, Gr){r,t,t'), is an even 
function of r, equation fll9l) implies that the transport coefficient Q;j/(r) vanishes. 

III. FORCED STOCHASTIC VELOCITY DYNAMICS 
A. Forced velocity dynamics for small Re 

We consider the simplest of dynamics for the velocity field by ignoring Lorentz forces, 
and assuming that the fiuid is stirred randomly by some external means. If the velocity 
fiuctuations have root-mean-squared (rms) amplitude v^-^s on some typical spatial scale i, 
the fiuid Reynolds number may be defined as Re = [vt-^s^/u), where u is the kinematic 
viscosity; note that Re has been defined with respect to the fiuctuation velocity field, not 
the background shear velocity field. In the limit of small Reynolds number (Re <^ 1), 
the nonlinear term in the Navier-Stokes equation may be ignored. Then the dynamics of 
the velocity field, v{X,t), is governed by the randomly forced, linearized Navier-Stokes 
equation. 





(26) 
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|: + "5^1^) ^ + "^^1^2 = -Vp + uV^v + / (27) 

f{X, t) is the random stirring force per unit mass which is assumed to be divergence-free 
with zero mean: V- / = and (/) = . The pressure variable, p, is determined by requiring 
that equation (1271) preserves the condition, V- v = 0. Then p satisfies the Poisson equation, 

V^P - (28) 

It should be noted that the linearity of the equations f l27p and (12 8 p implies that the velocity 
fluctuations have zero mean, (v) = . It is clear from equation (128|) that p is a non local 
function of the velocity field, so it is best to work in Fourier-space. Taking the spatial 
Fourier transform of equation (1271) . we can see that the Fourier transform of the velocity 
field, v{K,t), obeys. 



Or + "^-^j - - f ' = ^* '''' 



where fi{K, r) is the spatial Fourier transform of /j. It can be verified that the equation 
preserves the incompressibility condition K^Vm = 0. 

We can get rid of the inhomogeneous term, {K2d/dKi), in equation (12 9 p by transforming 
from the old variables {K,t) to new variables {k,t), through the Fourier-space shearing 
transformation of equation f l2T|) . First, we need to define new velocity and forcing variables, 
ai{k,t) and gi{k,t), respectively, by 

Vi{K,T) = G,{k,t,0)aiik,t) (30) 
/,(K,r) = G,{k,t,0)g^{k,t) (31) 
where Gu{k,t,0) is the Fourier-space viscous Green's function, defined by 



G^{k,t,t') = exp 



t 

2/ 



u / dsK'{k,s) 



(32) 



Noting the fact that K{k, s) = {ki — Ssk2, k2, k^) and K'^{k, s) = \K{k, s)p, the viscous 
Green's function can be calculated in explicit form as 
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Gy{k,t,t') = exp 



The Green's function possesses the following properties: 



G,{k,t,t') = G,{~k,t,t') 



G^{ki,k2,k3,t,t') = G^{ki,k2,-ks,t,t') 



Gu{k,t,t') = Gu{k,t, s) X Gu{k, s,t') , for any s 



Using the inverse transformation, 



(33) 



(34) 



Ki = ki - Stk2, K2 = ^2, = ks, 

and the fact that partial derivatives transform as, 



T = t 



(35) 



d 



d ^ ^ d 
oki 



dkj 



— = — + Sk — 

dr dt dki 



(36) 



equation ( 12^ leads to the following equation for the new velocity variables, ai{k,t): 



dtti [ K2Ki Si2 



(37) 



where K{k,t) = {ki — Stk2, k2, k^) and K'^{k,t) = |X(fc,t)p as given by equation (135|1 . 
It can be verified that equation fl37|l preserves the dot product, Kiai = 0. We also note 
that the dependence of the velocities, Vi{K,T) on the viscosity u arises solely through the 
Fourier-space Green's function. It is helpful to display in explicit form all three components 
of equation (157|1 : 



dai _ c^gf K1K2 



dt 



da2 
'dt 



das _ 2S ^ ^^^^ 



dt 



ai 



ai 



ai 



91 



92 



93 



(38) 
(39) 
(40) 
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Then equation ( 155]) can be solved to get an explicit expression for ai{k,t). When this is 
substituted in equations fl5^ and fHUl) . they can be integrated directly to obtain expressions 
for a2{k, t) and 03 (fc, t). The forced (or particular) solution, with initial condition ai{k, 0) = 
is 



= fdsg,{k,s)+ f ds[UKik,t)) - A,{Kik,s))]^^^^^g,{k,s) (41) 
Jo Jo ^± 



ai{ 



where Kj_ = K2 + = + kl = k"^, and the function, Aj, is defined as 



UK) = - ^ + ^ 



H2 



arctan 



(42) 



B. Velocity spectrum tensor expressed in terms of the forcing 

Our goal is to express the velocity spectrum tensor in terms of the statistical properties 
of the forcing. If the forcing is Galilean-invariant, then we must have. 



/,(K,r)/;(K',r')) = i2n)' 6ik - k')^,Uk,t,t') (43) 

where is the forcing spectrum tensor. We are now ready to use the dynamical solution of 
the last subsection. Using equations ( 130|) and (14T|) . Fourier-space, unequal-time, two-point 
velocity correlator is given by. 
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= G^{k,t,0)G^{k',t',0) [ ds [ ds' x 

Jo Jo 

xl^{g^{k,s)gUk',s')) + 

+ [A,{K{k,t)) - A,{K{k,s))]^^^^{g,{k,s)g*^^^^ 

+ [AUK{k',t')) - AUKik',s'))] ^'j^/^^ {g,{k,s)glik',s')) + 

+ [A,iK{k,t)) - A,iK{k,s))] [AmiK{k',t')) - AmiK{k',s'))] x 

K'^(k,s)K\k',s') , , ,,,, , 

K2^K''' ^1^^ ' ^ ^^^^ 

Using equations (13T|) and (H3|) . we write 



= ^- — ^ . ^ (f,{K{k,s),s)nK{k',s'),s'] 

^ (27r)^(5(fe-/e')$,™(fc,s,s') (45) 



G,(fc,s,0) G,ik',s',0) 
Using G^(fc,t,0)(G^(A;,s,0))-i = equations (glD, (05]) and ([22]) give, 
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Ujm{k,t,t') = [ ds [ ds' G^{k,t,s)G^{k,t',s') X 
Jo Jo 

x|$j^(fc,s,s') + 

+ [A,(K(fc,t)) - A,{Kik,s))] ^'l^^''U ,m{k,s,s') + 

+ [AUK{k,t')) - Am{K{k,s'))] ^'^^;''U ,,{k,s,s') + 

+ [A,{K{k,t)) - A,{K{k,s))] [AUK{k,t')) - AUK{k,s'))] x 
K^{k,s)K\k,s') 



<^n{k,s,s')} (46) 



When ^jm{k, t, t') is real, the forcing may be called non helical. Then equaton ( H6|) proves 
that the velocity spectrum tensor, Iljmik^t^t') is also a real quantity. In other words, non 
helical forcing of an incompressible fluid at low Re, in the absence of Lorentz forces, gives 
rise to a non helical velocity field. In this case, as we noted earlier, the velocity correlators 
Qjmi{f,t,t') and Cjmi{r,t,t') are odd functions of r and, Gj^{r,t^t') being an even function 
of r, equation flTIJl) implies that the transport coefficient, aji(r) vanishes. In other words, 
the a-effect is absent for non helical forcing at low Re and Rm, for arbitrary values of the 
shear parameter. This may not seem like a particularly surprising conclusion, but it is by 
no means an obvious one, because at high Re it may happen that Iljmik^t^t') is complex 
even when ^jm{k,t,t') is real. 

We now specialize to the case when the forcing is not only non helical, but isotropic and 
delta- correlated-in-time as well; in this case, 

$,„,(fc, s') = 6{s - s') P,m{K{k, s)) F (^^^) (47) 

where K{k,s) = \K{k,s)\, Kp = is the wavenumber at which the fluid is stirred, 
Pjm{K) = {6jm — KjKm/K'^) is a projection operator, and F{K/Kf) > is the forcing 
power spectrum. 

Substitute equation ( H7|) in ( H^ . and reduce the double-time integrals to a single-time 
integral using, 
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[ ds f ds'6{s-s')w{k,s,s') = f^dsw{k,s,s) (48) 
Jo Jo Jo 

where i< = Mm {t,t'). Then the velocity spectrum tensor, 



Ujm{k,t,t') = l'' dsG,{k,t,s)G,{k,t',s)F(^^^^] X 







X <^ Pjm{K{k, S)) + 



+ [Aj{K{k,t)) - A^{K{k,s))]^^^^P,m{K{k,s)) + 



K\k,s) 



+ [AUK{k,t')) - Am{Kik,s))] '^'^ P^,{K{k,s)) + 
+ [Aj{K{k,t)) - Aj{K{k,smArn{K{k,t')) - Am{K{k,s))] x 

x^^^Pn(K(/c,.))| (49) 

is completely determined when the forcing power spectrum, F{K/Kp), has been specified. 

Let an observer located at the origin of the laboratory frame correlate fluid velocities at 
time T = t and at time t' = t'. The two-point function that measures this quantity is given 

by, 

{Vj{0,T)Vm{0,T')) = Rjm{0,t,t') = J d^k Ujm{k,t,t') (50) 

It can be proved that, in the long time limit when t — )■ oo and t' — )■ oo, Rjm{0,t,t') 
is a function only of the time difference, {t — t'). The equal-time correlator, defined by 
Rjm{0,t,t) , is symmetric: Rjm{0,t,t) — Rmj{0,t,t). A related quantity is the root-mean- 
squared velocity, Vrms(^), defined by 

vLsit) = Rll{0,t,t) + R22{0,t,t) + i?33(0,t,t) (51) 

In the long-time limit, both Rj^{0,t,t) and Vrms(^) saturate due to the balance reached 
between forcing and viscous dissipation; let — limt_^oo 'f^rms(0- 
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We now define various dimensionless quantities: The fluid Reynolds number, Re = 
"^rms/ i^^p) ; the magnetic Reynolds number, Rm = f {rjKp) ; the Prandtl number, 
Pr = u /rj ; the dimensionless Shear parameter, Sh = S / {v'^^Kp) ■ 

For numerical computations, it is necessary to choose a form for the forcing power spec- 
trum. A quite common choice, used especially in numerical simulations, is forcing which is 
confined to a spherical shell of magnitude Kp. Therefore, whenever we need to choose a 
form for the forcing power spectrum, we take it to be. 



IV. PREDICTIONS AND COMPARISON WITH NUMERICAL EXPERIMENTS 

We have already established that the transport coefficient an = when the stirring is 
non helical. The other transport coefficient rjimi can be calculated by the following steps: 

(i) Computing the velocity spectrum tensor, Hjm, using equations fH9|l and fl52|) . 

(ii) Using this in equation to compute the velocity correlators Cjmi and -Djm- 

(iii) Substituting these correlators in the second of equations (fT9|) . 

We also seek to compare our analytical results with measurements of numerical simulations, 
which use the test-field method |5|. In this method, the mean-magnetic field is averaged 
over the coordinates Xi and X2. So we consider the case when the mean magnetic field, 
B = B{X^, r). The condition V- B = implies that is uniform in space, and it can be 
set to zero; hence we have B = {Bi, B2, 0). Thus, equation ( ITSjl for the mean EMF gives 
€ = {£1, £2, 0), with 

£^ = -V^J(r)J,; J = VXB = i-—^, —1, (53) 



. 9X3' 9X3' ^ 

where 2-indexed magnetic diffusivity tensor rjij has four components, {rjn, rju, 7721, ^722), 
which are defined in terms of the 3-indexed object rji^i by 

r]ij{T) = eij3r],ai{T); which implies that r]ii{T) = -r]i32ir), rii2{r) = rji^iir) (54) 
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Equation (1^ for S can now be substituted in equation ([T]). Then the mean-field induction 
becomes, 



The diagonal components, rju^r) and ?722(t), augment the microscopic resistivity, rj, whereas 
the off-diagonal components, ?7i2(t) and ri2i{T), lead to cross-coupling of Bi and i?2- 

A. The magnetic diffusivity tensor 

We now use our dynamical theory to calculate rjiji^r). From equations (IMl) and ffT^ . we 
have 



eij3(^ipm dr' j dV rs ^^(r, r, r) Cpmi{r, r, r') + 
+ 5,, ^ dr'^ dVG,(r,r,r') Z^33(r, r,r') 



(56) 



Thus the "D" terms contribute only to the diagonal components, rju and 7722- This is the 
expected behaviour of turbulent diffusion, which we now see is true for arbitrary shear. 
Using equation f l2^ . the velocity correlators Cpmi and D33 can now be written in terms of 
Hjm- After some lengthy calculations, the rjijlr) can be expressed in terms of the velocity 
spectrum tensor by. 



r/,,(r) = 2r,j dr' / d^fc r, r') (r - r') ^3 



6j2{ki - 5'r'A;2) - 6jik2 



6,1 {U23 - n32 - s{t - r')n3i} + 5,2 {n3i - ii,^} 



+ 



+ 6,, J dr' / d^fcG,(fc,r,r')n33 



(57) 



18 



K 

Eh 



CD 

K 



8 w 



N 
0) 

K 



0.3 
0.2 
0.1 



"T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r- 



1 — I — I — I — — I — I — I — I — — I — I — I — I — — I — I — I — I — — I — I — I — h 




FIG. 1: Plots of the saturated quantities rjt, r/^ and rj'^ for Re = Rm = 0.1 and Re = Rm = 0.5, 
corresponding to Pr = 1, versus the dimensionless parameter (— ShRe). The bold lines are for 
Re = Rm = 0.1, and the dashed lines are for Re = Rm = 0.5 . 



where 11/^ = ^im{k,T,T'), and the indices run over values 1 and 2. Here Gn{k,T,T') 
is the Fourier-space resistive Green's function defined in equation ([7]). The final step in 
computing rjiji^r) is to use equations ( H^ and for the velocity spectrum tensor. Him. 

The rjijlr) saturate at some constant values at late times; let us denote these constant 
values by 7]°° = rjijij — t- oo). If the mean magnetic field changes over times that are longer 
than the saturation time, we may use 7]"^ instead of the time-varying quantitites rjijij) 
in equation fl55|) . Looking for solutions B oc exp [Ar + i/^aXa], we obtain the dispersion 
relation, 



A 



2;iven m 



-1 ± 



riT 



+ e2 



(58) 



where the new constants are defined as. 
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12 3 4 5 

-Sj^Re 

FIG. 2: Plots of the saturated quantities rjt, rj'^ and t/^ for Re = 0.1 and Rm = 0.5, corresponding 
to Pr = 5, versus the dimensionless parameter (— ShRe). 



Vt = liVu +VT2), VT = V + Vt, e = ^{vTi - VT2) (59) 

Exponentially growing solutions for the mean magnetic field are obtained when the radicand 
in equation (l58i) is both positive and exceeds r]^ . 

From equations fl57j) . ([7]), fH9l) and fl52|) . it can be verified that the saturated values of 
the magnetic diffusivities, rj!^, have the following general functional form: 

„oo _ „ P 2 /ij(ShRe,Pr) 
- ^^^^ l + x(S.,Re,Pr) ' 

where the fij are dimensionless functions of two variables, and x is a dimensionless function 
of three variables. Figures (1-3) display plots of rjt, r]^ and rj^ , versus the dimensionless 
parameter (—ShRe). The scalings of the ordinates have been chosen for compatibility with 
the functional form displayed in equation (160|) above. These plots should be compared with 
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12 3 4 5 

-Sj^Re 

FIG. 3: Plots of the saturated quantities rjt, rj'^ and t/^ for Re = 0.5 and Rm = 0.1, corresponding 
to Pr = 0.2, versus the dimensionless parameter (— ShRe). 

Figure (3) of However, it should be noted that we operate in quite different parameter 
regimes; we are able to explore larger values of |Sh|, whereas jsj have done simulations for 
larger Re and Rm. The plots in Figure (la-c) are for Pr = 1, but for two sets of values of 
the Reynolds numbers; Re = Rm = 0.1, and Re = Rm = 0.5. Figure (2a-c) are for Re = 0.1 
and Rm = 0.5, corresponding to Pr = 5. Figure (3a-c) are for Re = 0.5 and Rm = 0.1, 
corresponding to Pr = 0.2. As may be seen from equation f l60p . the ratio, (Vu/v^)^ is a 
function only of the two dimensionless parameters, (ShRe) and Pr. In Figure (4) we plot 
this ratio versus (—ShRe) for all the cases considered in Figures (1-3). Some noteworthy 
properties follows: 

(i) We see that r]t is always positive. For a fixed value of (—ShRe), the quantity r]t/{r]TRe^) 
increases with Pr and, for a fixed value of Pr, it increases as (—ShRe) increases from 
zero (which is consistent with |[5|]), attains a maximum value near (—ShRe) ~ 2, and 
then decreases while always remaining positive. 
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FIG. 4: Plots of the ratio (t/i2/?/2i) versus the dimensionless parameter (— ShRe) for all the cases 
considered in Figures (1~3). The bold line is for the two cases corresponding to Pr = 1, the 
dashed-dotted line is for Pr = 5, and the dotted line is for Pr = 0.2. 

(ii) As expected, the behaviour of t]^ is more complicated. It is zero for (—ShRe) = 0, and 
becomes negative for not too large values of (—ShRe). After reaching a minimum value, 
it then becomes an increasing function of (—ShRe) and attains positive values for large 
(—ShRe). Thus the sign of is sensitive to the values of the control parameters. This 
may help reconcile, to some extent, the fact that different signs for 77^ are reported in 
12] and fl. 



iii) As may be seen, rj^ is always positive. This agrees with the result obtained in j^, |ll | 



and 



12| 



(iv) At first sight rj^ and t]^ appear to have quite different behaviour. However, closer 
inspection reveals certain systematics: as Pr increases, the overall range of values 
increases, while their shapes shift leftward to smaller values of (—ShRe). From equa- 
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tion (1^ . it is clear that the ratio {rj'^/r]^) is a function only of the two variables, 
(ShRe) and Pr. As Figure (4) shows, this ratio is nearly a linear function of (ShRe), 
whose slope increases with Pr. 

(v) The magnitude of the quantity, x(Sh5Re,Pr), that appears in equation flBUl) . is much 
smaller than unity. So rit/^riTRe^), r]^/ {riTRe^) and ri^/{riTRe^) can be thought of 
(approximately) as functions of (—ShRe) and Pr. This is the reason why, in Figure (1), 
the bold and dashed lines lie very nearly on top of each other. 



B. Implications for dynamo action and the shear current effect 

The mean magnetic field has a growing mode if the roots of equation (158|) have a positive 
real part. It is clear that the real part of A_ is always negative. So, for the growth of 
the mean magnetic field, the real part of A+ must be positive. Requiring this, we see from 
equation S5Eh that the condition for dynamo action is, 

J21 ^ IV2^ + — > 1. 61 

In Figure (5) we plot the last two terms, {VuV^ / Vt) ^-nd (e^/?7|.), as functions of (—ShRe), 
for all the four cases. Re = Rm = 0.1; Re = Rm = 0.5; Re = 0.1, Rm = 0.5 and 
Re = 0.5, Rm = 0.1. As may be seen, the magnitudes of both terms are much smaller 
than unity, so they are almost irrelevant for dynamo action. Hence, there is growth of the 
mean magnetic field only when the first term, [rj^S/rj'^K^) , exceeds unity. This is possible 
for small enough K^, so long as {f]^S) is positive. However, we see from Figures (1-3) 
that rj^ is always positive, implying that the product {i]^S) is always negative. Therefore 
the inequality of (!6T]) cannot be satisfied, and the mean-magnetic field always decays, a 



conclusion which is in agreement with those of [5|], [ll| and 12|. We can understand the 
above results more physically. Let us assume that {K^l is small enough, and keep only the 
most important terms in equation (155|) . Then we have. 



V21 + ■■■ , ^ - -^^i + • • • ' (62) 



dr dXl ' ' dr 



where we have used the saturated values of the magnetic diffusivity. If we now look for modes 
of the form B oc exp [Ar + iK^X^], we obtain the dispersion relation, A-t = ± K^^Jt]^^. So 
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FIG. 5: Plots of (e^/^/ji) and [rj^rj^/r]^^ versus the dimensionless parameter (— ShRe) for all the 
four cases considered in Figures (1-3). The bold lines are for Re = Rm = 0.1 ; the dashed lines are 
for Re = Rm = 0.5 ; the dashed-dotted lines are for Re = 0.1, Rm = 0.5 and the dotted lines are 
for Re = 0.5, Rm = 0.1. 



it is immediately obvious that A4. is real and positive — i.e. the mean magnetic field grows 
— only when the product (rj^S) is positive. However, this product happens to be negative, 
and the mean magnetic field is a decaying wave. 

The above results have direct bearing on the shear-current effect [lO^. This effect refers 
to an extra contribution to the mean EMF which is perpendicular to both the mean vorticity 
(of the background shear flow) and the mean current. From equation (l53i) . we see that, in 
our case, the relevant term is the contribution, —rj^Ji, to 82- As Figures (1-3) show, the 
diffusivity, 77^ is non zero only in the presence of shear, so the word shear refers to this. 
The word current refers to Ji, the cross-field component of the electric current associated 
with the mean-magnetic field 16|] . The shear-current effect would lead to the growth of the 
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mean magnetic field (for small enough K^), if only the product {r]^S) is positive. However, 
as we have demonstrated, this product is negative, so the shear-current effect cannot be 
responsible for dynamo action, at least for small Re and Rm, but for all values of the shear 
parameter. 



V. CONCLUSIONS 

Building on the formulation of [14], we have developed a theory of the shear dynamo 
problem for small magnetic and fluid Reynolds numbers, but for arbitrary values of the 
shear parameter. Our primary goal is to derive precise analytic results which can serve as 
benchmarks for comparisons with numerical simulations. A related goal is to resolve the 
controversy surrounding the nature of the shear-current effect, without treating the shear 
as a small parameter. We began with the expression for the Galilean-invariant mean EMF 



derived in [l^, and specialized to the case of a mean magnetic field that is slowly varying 
in time. This resulted in the simplification of the mean-field induction equation, from an 
integro-differential equation to a partial differential equation. This reduction is the first step 



to the later comparison with the numerical experiments of ^|. Explicit expressions for the 
transport coefficients, an and rjimi, were derived in terms of the two-point velocity correlators 
which, using results from [l4|, were then expressed in terms of the velocity spectrum tensor. 
Then we proved that, when the velocity field is non helical, the transport coefficient an 
vanishes; just like everything else in this paper, this result is non perturbative in the shear 
parameter. We then considered forced, stochastic dynamics for the incompressible velocity 
field at low Reynolds number. An exact, explicit solution for the velocity field was derived, 
and the velocity spectrum tensor was calculated in terms of the Galilean-invariant forcing 
statistics. For non helical forcing, the velocity field is also non helical and the transport 
coefficient an vanishes, as noted above. We then specialized to the case when the forcing is 
not only non helical, but isotropic and delta-correlated-in-time as well. We considered the 
case when the mean-field was a function only of the spatial coordinate and time r ; the 
purpose of this simplification was to facilitate comparison with the numerical experiments 
of |5|. Explicit expressions were derived for all four components, ?7ii(r), ri22{T) r]i2{T) and 
'721 (t), of the magnetic diffusivity tensor, in terms of the velocity spectrum tensor. Important 
properties of this fundamental object are as follows: 
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1. All the components of rjij are zero at r = 0, and saturate at finite values at late times, 



which we denote by rj. 



oo 



2. The off-diagonal components, 7712 and 7721, vanish when the microscopic resistivity 
vanishes. 

3. The sign of rj^ is sensitive to the values of the control parameters. This may help 



reconcile, to some extent, the fact that different signs for rj^ are reported in 12| and 

We derived the condition — the inequality fl6ip — required for the growth of the mean 
magnetic field: the sum of three terms must exceed unity. It was demonstrated that two 
of the terms are very small in magnitude, and hence dynamo action was controlled by the 
behaviour of one term. i.e. the mean magnetic field would grow if {rj^S/rj^K^) exceeds 
unity. This is possible for small enough K^, so long as {rj^S) is positive. However, we see 
from Figures (1-3) that 1]^ is always positive, implying that the product (rj^S) is always 
negative. Thus the mean-magnetic field always decays, a conclusion which is in agreement 



with those of 



0, [ni 



and 



12| . We then related the above conclusions to the shear-current 
effect, and demonstrated that the shear-current effect cannot be responsible for dynamo 
action, at least for small Re and Rm, but for all values of the shear parameter. In p], it is 
suggested that the dynamo action observed in their numerical experiments might be due to 
a fluctuating a-effect; addressing this issue is being the scope of our present calculations. 
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